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Abstract 

Foilless  diodes  used  to  produce  intense 
annular  relativistic  electron  beams  have  been 
simulated  using  the  time-dependent,  two-dimen¬ 
sional  particle-in-cell  code  C^UBE.  Current 
densities  exceeding  200  kA/cin  have  been 
obtained  in  the  simulations  for  a  5  MeV,  35  15 
diode.  Many  applications,  including  microwave 
generation,  collective  ion  acceleration  and 
high-density  plasma  heating  require  a  laminar 
electron  flow  in  the  beams.  The  simulation 
results  indicate  that  foilless  diodes  immersed 
in  a  strong  external  magnetic  field  can  achieve 
such  a  flow.  Diodes  using  technologically 
achievable  magnetic  field  strengths  (~100  kG) 
and  proper  electrode  shaping  appear  to  be  able 
to  produce  beams  with  an  angular  scatter  of 
less  than  35  mrad  at  the  current  densities  and 
energies  mentioned  above.  Scaling  of  the 
impedance  and  temperature  of  the  beam  as  a 
function  of  geometry,  magnetic  field  strength 
and  voltage  is  presented. 

Introduction 

Foilless  diodes  may  be  used  for  the  produc¬ 
tion  of  intense  annular  relativistic  electron  beams 
for  many  applications  including  microwave  genera¬ 
tion,  collective  ion  acceleration  and  high-density 
plasma  heating.1  Conventional  foil  diodes  have 
been  found  to  suffer  from  an  impedance  collapse 
when  plasma,  generated  by  electrons  striking  the 
anode  foil,  propagate  from  the  anode  to  the  cath¬ 
ode  thereby  electrically  shorting  the  diode.  This 
problem  is  eliminated  by  using  a  foilless  diode, 
thus  allowing  higher  current  densities  than  can  be 
obtained  with  a  foil  diode.  In  addition,  the  elec¬ 
tron  beam  generated  by  a  foilless  diode  is  not 

perturbed  by  passing  through  a  foil  nor  is  it  nec- 

2 

essary  to  replace  a  foil  for  repeated  operation. 

Although  there  has  been  some  investigation  of 
relativistic  electron  beam  generation  by  foilless 
diodes  a  firm  understanding  of  the  diode  has  been 


lacking.  We  have  analyzed  the  simple  diode  il¬ 
lustrated  in  Fig.  1  to  determine  the  scaling  of 
diode  impedance  and  beam  temperature  as  a  function 
of  geometry,  magnetic  field  strength  and  voltage. 
Some  investigators  have  assumed  that  the  foilless 
diode  impedance  is  determined  by  the  maximum  cur¬ 
rent  allowed  by  space  charge  in  the  drift  tube. 
Our  analysis  indicates  that  the  diode  impedance  is 
determined  by  the  equilibrium  that  the  beam  obtains 
which  is  not  necessarily  the  equilibrium  which 
gives  the  space-charge  limiting  current. 

Impedance  Model 

Because  most  applications  require  a  beam  with 
laminar  flow  it  is  useful  to  model  the  beam  formed 
by  the  foilless  diode  by  the  cold  fluid  equations . 
In  an  azimuthally  symmetric,  axially  homogeneous 
equilibrium,  the  equations  describing  the  beam 


depend  only  on  the  radial  coordinate,  r. 

tions  to  be  solved  are 

The  equa 

mc2/e  y|3e2/r  =  -  PzB0 

(1) 

dBz/dr  =  4nne(30 

(2) 

d(rB0)/dr  =  -4nnerPz 

(3) 
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d(rEr)/dr  =  -4nner  (4) 

where  m  and  e  are  the  mass  and  charge  of  the  elec¬ 
tron.  The  only  nonzero  fluid  variables  are  the 
density  n  and  the  axial  and  azimuthal  fluid  veloc¬ 
ities  p  and  Pe  (divided  by  the  speed  of  light  c) . 
The  nonzero  fluid  variables  are  the  radial  elec¬ 
tric  field  E^.,  the  azimuthal  magnetic  field  Bg,  and 
the  axial  magnetic  field  The  relativistic 

factor  is  denoted  by  y.  Because  the  cathode  is  an 
equipotential  surface,  conservation  of  energy 
assumes  the  following  form: 

dy/dr  =  -eEr/mc2  (5) 


Because  there  are  only  five  equations  and  six 
unknowns  another  condition  must  be  specified.  A 
condition  which  leads  to  an  analytically  tractable 
solution  of  the  equilibrium  equations  and  which 

becomes  increasingly  better  satisfied  at  larger 

g 

energies,  is  to  choose  |5z  to  be  independent  of  r. 

2  -Is 

Defining  the  following  quantities  Yj |  =  (1  "  Pz  )  * 
and  y,  =  Y/Yj|  an  equation  for  y^  can  be  found 
whose  solution  is  given  in  terms  of  elliptic  Jacobi 

9 

functions.  The  total  beam  current  v,  measured  in 

3 

units  of  me  /e  is  given  by  this  model  as 

V  =  [(Y||  -  l)(yf0  -  1)(Y*0  +  a2)]%/2  (6) 


where  y^Q  is  y^  evaluated  at  the  outside  edge  of 
the  beam,Ro,  and  a  is  an  arbitrary  constant.  De¬ 
fining  uj  =  eB  (R  )/mc2  we  find 

CO  2  0 


Vco/c  =  (YIo  '  1)%  +  *io(vlo  +  a2)%  (7) 

and  In  R  /Ri  =  F((|>,k)/(a2  +  l)*5  W 


where  R.  is  the  inside  beam  radius,  $  =  Cos  *(y  7) 
and  k  =  a/ (a^  +  l)"1  and  F(4>,k)  is  the  incomplete 
elliptic  integral  of  the  second  kind. 

In  addition  to  Eqs.  (6)-(8),  we  require  that 
the  total  energy  of  the  electrons,  kinetic  and 
potential,  be  equal  to  the  potential  drop  between 
the  anode  and  cathode.  Thus, 


y  =  y, .  y,  +  2u/p  In  R  /R  (9) 

*a  *  I  I  io  z  a  o 

where  R  is  the  anode  radius  (see  Fig.  1).  The 

3 

relativistic  factor  that  the  electrons  would  have 


upon  reaching  the  anode  is  denoted  by  y  .  Voronin, 
et  al.  have  used  these  equations  and  additional 
assumptions  to  find  the  space-charge  limited  cur- 

g 

rent  as  a  function  of  magnetic  field.  However, 
there  is  no  a  priori  reason  to  assume  that  the 
beam  produced  by  the  diode  will  be  launched  into  an 
equilibrium  which  will  transmit  the  maximum  current. 

If  the  applied  external  magnetic  field  pene¬ 
trates  the  cathode  then  conservation  of  canonical 
angular  momentum  takes  the  form: 

(y2  -1)*  =  (R  iu  /c  -  R2  u>  /R  c)/2  (10) 

Mo  o  co  c  c  o 

2 

where  tuc  =  eBg/cm  and  Bg  is  the  applied  magnetic 

field.  The  cathode  radius  is  denoted  bv  R  .  If 

'  c 

in  addition  we  assume  that  the  laminar  electron 
flow  is  along  the  self-consistent  magnetic  field 
lines,  then  the  flux  between  the  axis  and  the  outer 
edge  of  the  beam  is  equal  to  the  applied  flux  be¬ 
tween  the  axis  and  the  cathode  radius.  On  the  time 
scale  of  most  experiments,  the  magnetic  field  pro¬ 
duced  by  the  beam  cannot  diffuse  through  the  anede 
wall.  Therefore,  the  flux  between  the  outer  edge 
of  the  beam  and  the  anode  will  be  equal  to  the 
applied  flux  between  the  cathode  and  anode.  These 
conditions  may  be  written  as 

R  2U)  / c  =  2  R  (y  2  -  \)h  +  R.(l  +  a2)%  (11) 

C  C  O  1 0  1 

and  tn  (R2  -  R2)  =  tu  (R2  -  R2)  .  (12) 

co  a  o  c  a  c 

Equations  (6)-(12)  form  a  complete  set  of  equa¬ 
tions  which  can  be  solved  (numerically)  to  determine 
the  impedance  of  the  foilless  diode.  In  order  to 
insure  laminar  flow,  it  is  necessary  to  apply  a 
large  external  magnetic  field.  Therefore  a  useful 
approximation  can  be  obtained  by  taking  the  infinite 
magnetic  field  limit.  One  then  finds  that  the  beam 
becomes  infinitesimally  thin  with  radius  R^  and  that 
the  beam  approaches  a  nonrotating  equilibrium.  The 
diode  impedance  in  this  limit  becomes 

Z  =  15(ya-l){[ya/(l+4  In  Ra/Rc) ]2-l}-%  Cl.  (13) 

should  be  noted  that  this  formula  is  invalid 
for  low  voltage,  probably  owing  to  our  assumption 
of  being  independent  of  r. 
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Diode  Simulations 

A  two-dimensional  relativistic  time-dependent 
particle-in-cell  simulation  code,  CCUBE,  has  been 
used  to  test  the  impedance  model  and  gain  insight 
into  the  parameters  affecting  beam  quality.^  An 
emission  algorithm  in  the  code  emits  charge  from 
the  cathode  surface  at  a  sufficient  rate  to  satisfy 
the  space-charge  limited  emission  boundary  condi¬ 
tion,  i.e.,  the  electric  field  normal  at  the  cathode 
is  zero.  The  diode  simulations  were  run  with  a 
transverse  electromagnetic  (TEM)  wave  launched  from 
the  left  in  Fig.  1  onto  the  coaxial  transmission 
line.  By  not  allowing  the  first  few  cells  to  emit, 
one  can  control  the  impedance  of  the  driver  to  the 
diode,  which  in  all  cases  was  taken  to  be  37  Q. 
Typically  the  length  of  the  simulation  region,  L, 
was  5  to  10  cm.  Impedances  and  beam  parameters 
are  measured  when  the  system  consisting  of  the 
transmission  line  driver  with  the  diode  load  had 
reached  steady  state.  At  this  time  the  voltage 
on  the  diode,  V,  is  given  by 

V  =  2VwZ/(Zq  +  Z)  (14) 

where  Z  is  the  diode  impedance,  ZQ  is  the  transmis¬ 
sion  line  impedance  and  is  the  voltage  of  the 
TEM  wave  launched  onto  the  line.  Diagnostics  in 
the  code  include  voltage  and  particle  current 
probes,  Rogowskii  Coils  as  well  as  impedance  probes 
located  at  several  axial  positions.  At  the  end  of 
the  simulation  region  diagnostics  include  Faraday 
Cups,  calorimeters  and  density,  mean  velocity,  and 
temperature  measurements  as  a  function  of  radial 
position. 

Simulation  Results 

From  Eq.  (13)  we  see  that  for  large  applied 
magnetic  fields  the  diode  impedance  depends  only 
on  the  voltage  and  the  ratio  of  the  anode  to  cathode 
radius.  Figure  2  shows  the  results  of  several 
simulations  performed  with  a  Vw  =  5 . 1  MV  and  a 
cathode  radius,  Rc  =  1  cm.  Because  of  the  impedance 
mismatch,  the  voltage  across  the  diode  varied  from 
3.5  to  6.0  MV  in  accordance  with  Eq.  (14).  The  open 
circles  represent  simulations  performed  with  an 
applied  magnetic  field  of  100  kG  and  an  A-K  gap, 

5,  (see  Fig.  1)  of  0.4  cm.  The  triangle  represents 
a  run  with  the  same  parameters  but  with  5  =  0.2  cm. 


Fig.  2.  Current  versus  ratio  of  anode  to  cath¬ 
ode  radius  for  various  foilless  diodes . 

The  dashed  line  is  from  E'q.  (13).  The 
solid  line  is  the  space-charge  limit. 

Two  runs  at  40  kG  and  6  =  0.4  cm  are  denoted  by 
inverted  triangles.  The  square  designates  a  run  at 
55  kG  in  which  the  anode  wall  is  continued  straight 
at  the  original  transmission  line  radius,  R^,  of 
1.85  cm  so  that  5  ■*  The  dashed  line  is  obtained 
from  Eq.  (13).  The  solid  line  is  the  space  limiting 

current  for  the  infinitesimally  thin  beam  in  the 

14 

infinite  magnetic  field  limit.  The  simulation 
data  in  all  cases  lies  well  below  the  space-charge 
limiting  current  and  rather  close  to  the  current 
given  by  the  impedance  formula  of  Eq.  (13).  All 
the  simulations  were  performed  with  a  cathode  tip 
thickness,  s  in  Fig.  1,  of  0.135  cm  except  the  run 
at  55  kG  which  had  e  equal  to  R^.  At  100  kG  the 
beam  thickness  was  found  to  be  approximately 
0.03  cm,  thus  it  is  unlikely  that  much  effect  would 
be  found  for  e's  larger  than  this  value.  These  very 
thin  beams  can  yield  current  densities  exceeding 
700  kA/cm2. 

Because  Eq.  (13)  was  obtained  for  the  infinite 
magnetic  field  limit,  it  is  desirable  to  determine 
the  effects  of  finite  magnetic  field.  Figure  3 
shows  the  results  of  a  series  of  simulations  per¬ 
formed  with  V  =  5.1  MV,  R  =1.27  cm,  6=0. 4  cm, 
w  a 


7; 


and  £=0.135  cm.  The  dashed  curve  is  obtained  from 
the  numerical  solution  of  Eqs .  (6)-(12).  The  solid 
curve  is  the  space-charge  limited  diode  theory  of 
Voronin,  et  al.10  The  diode  operates  well  below 
the  space-charge  limit  and  again  agrees  well  with 
the  laminar  flow  impedance  model.  The  beam  strikes 
the  anode  wall  at  25  kG  and  for  all  values  of 
applied  field  below  this  value  the  transmitted  elec¬ 
tron  current  gradually  diminishes  owing  to  current 
loss  to  the  anode.  Although  the  diode  impedance 
does  not  vary  much  with  magnetic  field,  the  temper¬ 
ature  as  measured  by  the  mean  annular  scatter  of 
electrons  around  the  beam  propagation  direction  was 
found  from  the  simulations  to  vary  from  200  mrad  at 
27  kG  to  less  than  60  mrad  at  100  kG.  The  large 
magnetic  field  makes  it  more  difficult  for  the  elec¬ 
trons  to  cross  field  lines  and  create  temperature 
by  mixing. 

One  would  expect  that  as  6  is  increased  to 
larger  values  that  the  beam  produced  by  the  foil¬ 
less  diode  would  come  to  equilibrium  before  it 

"sees"  the  reduced  anode  radius  E  .  The  existence 

a 

of  this  effect  is  verified  by  the  series  of  simula¬ 
tions  shovra  in  Fig.  4.  The  parameters  of  the 
simulations  include  BQ=100  kG,  Vw=5.1  MV,  Ra=1.5  cm, 
and  £=0.135  cm.  The  upper  dashed  curve  was  calcu¬ 
lated  from  Eq.  (13).  The  lower  dashed  curve  was 
also  calculated  from  Eq.  (13)  but  with  Ra  equal  to 

outer  radius  of  the  transmission  line  feed  R,.  As 

d 

seen  from  the  data,  the  diode  impedance  makes  a 
sudden  transition  from  an  equilibrium  with  an  anode 
radius  of  Ra  for  small  6  to  an  equilibrium  with 
anode  radius  at  large  6.  The  value  of  6  at  the 
transition  point  is  roughtly  one-half  of  the  cathode 
radius  for  this  case.  The  actual  transition  point 
is  probably  governed  by  the  distance  the  beam  must 


Fig.  3.  The  effect  of  finite  magnetic  field 
on  diode  impedance. 


Fig.  4.  Current  versus  A-K  gap,  £,  for  the 
foilless  diode. 

propagate  from  the  cathode  tip  to  reach  equilibrium 
and  must  certainly  depend  upon  the  current  density. 

For  large  values  of  6,  the  beam  gains  kinetic 
energy  as  it  approaches  the  region  in  which  the 
anode  radius  is  reduced  to  making  it  stiffer  and 
less  likely  to  phase  mix.  Simulations  have  shown 
beam  temperatures  below  25  mrad  for  this  scheme, 
which  is  near  the  numerical  resolution  of  the  code. 
Use  of  these  ideas  in  diode  design  show  promise  for 
producing  very  laminar  beams . 
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